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ABSTRACT 

A model with non-uniform flow in the longitudinal direction is pro- 
posed for the relativistic heavy-ion collisions and compared with the 14.6 
A GeV/c Si-Al and 10.8 A GeV/c Au-Au collision data. The stronger 
influence of transparency on the distribution of heavier produced particles 
and the larger stopping in the heavier collision system are accounted for 
by using a new geometrical parameterization picture. The central dips 
in the proton and deuteron rapidity distributions for Si-Al collision are 
reproduced. 
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I. Introduction 

The experimental finding that colliding nuclei are not transparent but undergo 
a violent reaction in central collisions represents one of the major motivations for 
the study of ultra-relativistic heavy ion collisions at the CERN/SPS, BNL/AGS and 
also at the future BNL/RHIC and CERN/LHC. Of central importance is the ability of 
understanding to what extent the nuclear matter has been compressed and heated. 

The use of thermal models to interpret data on particle distribution from nuclear 
collisions is motivated by the hope and expectation that in collision between sufficiently 
large nuclei at sufficiently high energies a state of exicted nuclear matter close to local 
thermodynamic equilibrium can be formed, allowing us to study the thermodynamics 
of QCD and the possible phase transition from a hadronic gas to a QGP. Since we do 
not yet reliably know how big the collision system and how large the beam energy have 
to be for this to occur (if at all), thermal models should be considered as an importent 
phenomenological tool to test for such a behavior. 

The study of collective flow in high energy nuclear collisions has attracted increasing 
attentions from both experimental' 1 ' and theoretical' 2 ' point of view. The rich physics 
of longitudinal and transverse flows is due to their sensitivity to the system evolution 
at early time. The expansion and cooling of the heated and highly compressed matter 
could lead to a considerable collectivity in the final state. Due to the high pressure, 
particles might be boosted in the transverse and longitudinal directions. The collective 
expansion of the system created during a heavy-ion collision implies space-momentum 
correlation in particle distributions at freeze-out. 

The experimental data on the rapidity distributions of produced particles in 14.6 
A GeV/c Si-Al collisions has been ultilized to study the collective expansion using a 
cylindrical-symmetric longitudinal flow model (CSFM)' 3 ' 4 '. In this model, the distri- 
bution of final-state particles comes from the superposition of a number of fire-balls, 
distributed uniformly within the kinematical limit along the longitudinal (rapidity y) 
axis. The model-results fit well the experimental distribution of pion, but are too nar- 
row in the case of heavier particles, such as proton and deuteron. In particular, the 
central dip, which can be clearly seen in the distribution of proton, is not reproduced. 

More recently, E877 Collaboration' 5 ' has published their data for 10.8 A GeV/c Au- 
Au collisions, which provides a good chance to compare the stopping power in collision 
systems of different sizes. A possible central peak of the rapidity distribution of proton 
at around mid-rapidity, which was obtained through extrapolating the experimental 
data to mid-rapidity using RQMD model' 6 ', has been taken as an evidence for the in- 
creasing of stopping power, but the reliability of this extrapolation is model-dependent. 

Stopping and transparency describe the same physical aspect in relativistic heavy 



ion collision from two opposite sides. In order to settle a reliable model for relativistic 
heavy ion collision, this aspect must be taken into account carefully. The assumption 
of uniformly distributed fire-balls in the CSFM model is a crude one. It does not 
account for the memory of the fire-balls on the motion of the incident nuclei. In the 
present paper we propose a non-uniform longitudinal flow model (NUFM) using a new 
geometrical parameterization picture^ to describe the nuclear stopping/transparency 
in a more proper way. The stronger influence of transparency on the distribution of 
heavier produced particles (proton and deuteron) as well as the larger stopping in the 
heavier collision system (Au-Au) are described by using this picture. The central dips 
in the proton and deuteron rapidity distributions for Si-Al collisions are reproduced, 
and at the same time the central peak in the pion rapidity distributions is maintained. 

In section II the non-uniform longitudinal flow model (NUFM) with a geometrical 
parameterization picture is formulated. The results of the model are given and com- 
pared with the experimental data in section III. A short summary and conclusion is 
given in section IV. In order to avoid the complexity associated with the production 
of strange particles and concentrate on the expansion of the system, we will discuss in 
this paper only normal non-strange particles - - pions, protons and deuterons. 



II. Non-uniform longitudinal flow 



Firstly, let us briefly recall the fireball scenario of relativistic heavy ion collisions. 
Since the temperature at freeze-out exceeds 100 MeV, the Boltzmann approximation 
is used. Transformed into rapidity y and transverse momentum p t this implies^: 

E ^ oc Ee(~ E ^ = m t coshy e (-^^ y /T) (1) 
d 6 p 



Here m t = y m? + pi is the transverse mass, m is the mass of the produced particles 
at freeze-out. 

The rapidity is defined as y = tanh" 1 (pi / E), where p\ is the longitudinal momentum 
of the produced particle. Substituting into Eqn.(l) and integrating over m t , we get the 
rapidity distribution of the isotropic thermal source, 

^«ga +2i+1 aw ( 2 ) 

where £ — T/m cosh y. 

Eqn's.l and 2 give the isotropic thermal distribution. As mentioned in Ref. [4], 
for pions with m ~ T, this distribution is close to that of massless particles, i.e. 



proportional to 1/ cosh 2 y. For heavier particles isotropic emission implies a strongly 
narrow distribution, in contradiction to the experimental findings, see e.g. the dashed 
lines in Fig. 6. 

The measured momentum distribution of the final-state particles is certainly aniso- 
tropic. It is privileged in the direction of the incident nuclei. This is because the 
produced hadrons still carry their parent's kinematic information, making the longi- 
tudinal direction more populated than the transverse ones. The simplest way' 3,4 ^ to 
account for this anisotropy is to add up the contributions from a set of fire-balls with 
centers located uniformly in the rapidity region [—yo,yo], as sketched schematically in 
Fig.l. The corresponding rapidity distribution is obtained through changing the £o in 
Eqn.(2) into £ = T/m cosh(y — y') and integrating over y' from — y to y&. 

where £ = T/m cosh(y — y'). Equivalently, we can also use the angular variable 6 
defined by = 2 tan^ 1 exp(— y'), and change the integration variable in Eqn.(3) to 6, 

^ = eKm > T /*- J» (, + U + . g ) 
ay Je min sinH y mcosh(y — y') m 2 cosh(y — y')J 

x exp(-m cosh(y - y') /T) , (4) 

cf. the solid circle and lines in Fig. 2. 

This simple approach fits the rapidity distribution of pions well but fails to re- 
produce the central dip in heavier produced particles, which is clearly seen in the 
experimental distributions of protons and deuterons. 

Note that in this CSFM model the distribution of fire-balls in the longitudinal direc- 
tion of phase space is uniform which does not account for the interaction of the incident 
nuclei when they pass through each other properly. This is a crude approximation. A 
more reasonable assumption is that the fire-balls keep some memory on the motion of 
the incident nuclei, and therefore the distribution of fire-balls, instead of being uniform 
in the longitudinal direction, is more concentrated in the direction of motion of the 
incident nuclei, i.e. more dense at large absolute value of rapidity. A parametrization 
for such a distribution can be obtained by using an ellipse like picture on emission 
angle distribution, as shown in Fig. 3. In this scenario, the emission angle is 

9 = tan -1 (e tan 9), (5) 

where the parameter e(0 < e < 1) represents the ellipticity of the introduced ellipse 
which describes the non-uniformity of fire-ball distribution in the longitudinal direction, 
as sketched in Fig. 4. 



Subsituting Eqn.(5) together with y e = — In tan(0/2) into Eqn.(3), a distribution 
function Q{6) for the emission angle of non- uniform flow is introduced in the integration 

= «^ f , + h 2 ( r , + , 2 X 

dy ie min sine/ \ mcosh(y-y c ) m 2 cosh (y - y e ) / 

xexp(-mcosh(y-y c )/T), (6) 

y. = - In tan(0/2) , Q(9) = \ -. (7) 

Ve + tarn \ cos#| 

Here 6 min = 2tan _1 (e _ye0 ), # max = 2 tan" 1 (e J/c0 ). y c0 is the rapidity limit which confines 
the rapidity interval of longitudinal flow. 

Changing the integration variable in Eqn.(6) back to y e , the rapidity distribution 
can be rewritten as follows: 



d "— = eKm 2 T jT fl + + ^ 



dy i-2/eo \ mcosh(y-y c ) m 2 cosh (y - y c ) / 

x exp(-m cosh(y - y c )/T)p(y e )dy e . (8) 



1 + sinh 2 (y e 



(9) 



\ 1 + e 2 sinh 2 (y c ) 

Here p(y e ) is the distribution function of non-uniform flow in the longitudinal direction. 
It can be seen from fig. 5, that the larger is the parameter e, the flatter is the distribution 
p(y e ) and the more uniform is the longitudinal- flow distribution. When e — > 1, the 
longitudinal-flow distribution is completely uniform (p(y e ) — > 1), returning back to the 
CSFM model. 



III. Comparison with experiments 

We now proceed to compare the model-results with the AGS data. 

The rapidity distributions of pion, proton and deuteron for 14.6 A GeV/c Si-Al 
collisions^ 9 ' 10 ' are given in Fig.6 (a, b and c). The dashed, dotted and solid lines (band) 
correspond to the results from isotropical thermal model, uniform longitudinal flow 
model (CSFM) and the non-uniform longitudinal flow model (NUFM) respectively. The 
rapidity limit y c o and the ellipticity e used in the calculation are listed in Table I and 
illustrated in Fig. 7. The rapidity limit y used in the CSFM model of Ref. [4] is also 
listed in Table I for comparison. The parameter T is chosen to be 0.12 GeV following 
Ref. [4]. 

Note that the distribution of the light particle (pion) is insensitive to the ellipticity 
e. Changing e from 0.35 to 0.7 results in a narrow band shown in Fig. 6(a). Since 



pions are produced through the interaction of colliding nuclei, they have less memory 
on the motion of the incident nuclei before interaction. Therefore, physically the value 
of ellipticity e for pion should be bigger than that for proton and deuteron, but the 
exact value cannot be fixed through fitting the model results to the experimental data. 



Table I The value of model-parameters 



Parameter 


Si-Al Collisions 


Au-Au Collisions 


71 


P 


d 


71 


P 


e 


0.35 - 0.7 


0.52 


0.56 


0.35 - 0.7 


0.58 




1.35 


1.35 


1.35 


1.05 


1.05 




1.15 


1.15 


1.15 







It can be seen from the figures that the NUFM model reproduces the central dip 
of the rapidity distribution of heavier particles (proton and deuteron) in agreement 
with the experimental findings, while for light particles (pions) there is a peak instead 
of dip at central rapidity. Note that the appearance or disappearance of central dip 
is insensitive to the rapidity limit y e o but depends strongly on the magnitude of the 
ellipticity e and the mass m of the produced particles. For the heavier particles (proton 
and deuteron) a central dip appears for e < 0.8, but for light particles (pions) there is 
no dip even when e is as small as 0.35, cf. Table I and Fig. 6. 

On the other hand, the width of the rapidity distributions is mainly controlled by 
the amplitude y e o of the longitudinal flow. A single value (1.35) of y c o can account for 
the wide distribution of heavier particles (protons and deuterons) and at the same time 
fit the pion-distribution well. The difference in the width of the dN/dy distribution 
originates mainly from the difference in the mass of the particles' 3 ' 4 '. 

In Fig. 8 are shown the rapidity distributions of pions and protons for Au+Au 
collisions at 10.8 A GeV/c^. The dotted and solid lines (band) correspond to the 
results of CSFM and NUFM models (with parameters listed in Table I) respectively. In 
the calculation the CSFM results we have also used the NUFM model with the same 
rapidity limit y C Q listed in Table I but with ellipticity e = 1. The histogram is the 
result from the RQMD model. 

It can be seen from Fig. 8 that in the NUFM model there is a shallow dip in the 
central rapidity of the distribution of proton, instead of a central peak as predicted 
by the CSFM model. However, the presently available experimental data are restricted 
to the large rapidity. The peak at central rapidity is the extrapolation of data using 
RQMD and is model dependent. It is interesting to see whether the prediction of a 
central dip (plateau) or a central peak will be observed in future experiments. 

Comparing the parameter values for Si-Al (smaller colliding nuclei) and Au-Au 
(larger colliding nuclei) collisions listed in Table I, it can be seen that the rapidity limit 



y c0 is smaller and the ellipticity e for proton is bigger for the larger colliding nuclei than 
for the smaller ones. Both of these two show that the hadronic system formed from 
the larger colliding nuclei is less spread out in rapidity, i.e. there is stronger nuclear 
stopping in the collision of larger nuclei. 

IV. Summary and Conclusions 

In high energy heavy-ion collisions, due to the transparency of nucleus the partic- 
ipants will not lose the historical memory and the produced hadrons will carry some 
of their parent's memory of motion, leading to the unequivalence in longitudinal and 
transverse directions. So it is reasonable to assume that the flow of produced par- 
ticles is privileged in the longitudinal direction. This picture has been used by lots 
of models' 3,8 ]. Here we would mention two thermal and hydrodynamic models, one is 
the boost-invariant longitudinal expansion model postulated by Bjorken ' 8 1 which can 
explain such an anisotropy already at the level of particle production in hadron-hadron 
collisions. This model has been formulated for asymptotically high energies, where the 
rapidity distribution of produced particles establishes a plateau at midrapidity. The 
second model is the CSFM postulated first by Schnedermann, Sollfrank and Heinz® 
which accounts for the anisotropy of longitudinal and transverse direction by adding 
the contribution from a set of fire-balls with centers located uniformly in the rapidity 
region [-yo,Vo\ i n the longitudinal direction, sketched schematically in Fig's. 1 and 2 as 
dashed circles. It can account for the wider rapidity distribution when comparing to 
the prediction of the pure thermal isotropic model but fails to reproduce the central 
dip in the proton and deuteron rapidity distributions. 

In this paper, we argue that the transparency/stopping of relativistic heavy ion 
collisions should be taken into account more carefully. It will not only lead to the 
anisotropy in longitudinal-transverse directions, but also render the fire-balls (especially 
for those of proton and deuteron) concentrate more in the direction of motion of the 
incident nuclei. A non-uniform longitudinal flow model is proposed, which assumes 
that the centers of fire-balls are distributed non-uniformly in the longitudinal phase 
space. A parameter e is introduced through a geometrical parametrization which can 
express the non-uniformity of flow in the longitudinal direction, i.e. the centers of fire- 
balls of proton and deuteron prefer to accumulate in the two extreme rapidity regions 
(y e ~ ±2/ c0 ) in the c.m.s. frame of relativistic heavy-ion collisions, and accordingly the 
distribution is diluted in the central rapidity region (y c 0). 

It is found that the depth of the central dip depends on the magnitude of the 
parameter e and the mass of produced particles, i.e. the non-uniformity of longitudinal 
flow which is described by the parameter e determines the depth of the central dip for 



heavier particles. On the other hand, the central peak in the pion distribution turns 
out to be insensitive to the value of e and can be well reproduced from this model 
simultaneously with the central dip in the proton and deuteron rapidity distributions. 

Through comparing the feature of collision systems of different size, it is found that 
the maximum flow velocities are smaller for the heavier collision system than for the 
lighter ones, which suggests, together with the larger e, a larger stopping in the bigger 
collision system. 
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Figure captions 



Fig.l Schematic sketch of the distribution of fire-balls in the uniform flow model 
(CSFM). 

Fig. 2 Schematic sketch of the emission angle 6 (solid circle and lines) and the corre- 
sponding distribution of fire-balls (dashed circles) in the uniform flow model (CSFM). 

Fig. 3 Schematic sketch of the emission angle 9 in the non-uniform flow model (NUFM). 

Fig. 4 Schematic sketch of the distribution of fire-balls in the non-uniform flow model 
(NUFM). 

Fig. 5 The distributions p(y e ) of the center of fire-balls as a function of y e . When 
e — > 1, the non- uniform distribution turns to the uniform distribution p(y e ) — > 1. 

Fig. 6 Rapidity distributions for central 14.6 A GeV/c Si+Al collisions. Open circles 
are the experimental data for Si+Al collisions taken from Ref.'s [9,10]. Dashed, dotted 
and solid lines are the distributions from the isotropical thermal model, cylindrical- 
symmetric longitudinal flow model (CSFM) and non-uniform longitudinal flow model 
(NUFM) respectively. Fig's, (a), (b) and (c) are the pion, proton and deuteron distribu- 
tions respectively. The ellipticity e for pion is within a range [0.35, 0.7] and therefore 
the NUFM results are present as a band. The temperature T = 0.12 GeV. 

Fig. 7 The fire-ball distribution functions p(y e ) verus rapidity y e in the non-uniform 
flow of different particles for the collisions of Si+Al at 14.6 A GeV/c (a) and Au+Au 
at 10.8 A GeV/c (b) respectively. For pion the distribution functions p(y e ) are only 
ploted for the two boundaries of the used region 0.35 < e < 0.7. 

Fig. 8 Rapidity distributions for pions (a) and protons (b) in central Au+Au collisions 
at 10.8 A GeV/c. Full circles represent the experimental data taken from Ref. [5]. Open 
circles are the reflection of the data. The solid line is our calculation using the NUFM 
model. The ellipticity e for pion is within a range [0.35, 0.7] and therefore the results 
for pion are present as a band. The histogram shows the results from RQMD and 
the dotted lines are the results from the prediction of CSFM model. The temperature 
T = 0.14 GeV. 
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Fig. 7 
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